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In this work we study the anomalous top quark-gluon couplings Chromoelectric Dipole Moment 
(CEDM) and Chromomagnetic Dipole Moment (CMDM) in a general THDM with CP violation. We 
hnd that this model provides an important contribution from the Yu Yukawa coupling that needs to 
be taken into account. The prediction for CMDM and CEDM obtained are —0.03 < Akt < —0.005 
and |Adt| < 0.005, respectively. 
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I. INTRODUCTION 

The recent discovery of a Higgs Boson with a mass of 125 GeV in the Large Hadron Gollider (LHC) [U [2] has 
confirmed that the Standard Model (SM) is the theory that correctly describes electroweak interactions. However, 
in the context of the SM there are many unsolved problems, among them one can mention: (a) the fermion mass 
hierarchy, where the top quark is much heavier than other fermions, being the heaviest particle with a mass around 
the symmetry breaking scale, (b) CP symmetry breaking in the CKM matrix cannot explain the matter antimatter 
problem and the fermionic electric moments.This issues make very important the study of new sources of CP violation 
beyond the SM. 

On the other hand, top quark physics is a relevant scenario for the study of physics beyond the SM 0. In the 
next LHC run, at ^/s = 14 TeV, millions of top pairs will be produced, giving a great opportunity for the study of 
top quark properties , including its couplings to the gauge fields and providing an excellent scenario for new physics 
searches. 

In the SM, top quark magnetic and chromomagnetic dipole moments (MDM and CMDM) are induced at one 
loop level; on the other hand, top quark anomalous electric and chromoelectric dipole moments (EDM and CEDM) 
appear only at 3 loop level, arising from the complex phase in the CKM matrix. This couplings are important in 
the study of new physics; in order to generate them, extended models with new sources of CP violation are required. 
This anomalous couplings sensitively affect top pair production in pp collisions. Indirect bounds to the top quark 
MDM have been found from the bottom quark radiative decay 6 —>■ sq HIS], and also semileptonic B meson decays 
{B —>■ Kl'^l~) have been used to improve indirect bounds to MDM and EDM [5]. 

Since anomalous top quark couplings (ttg) affect the top production, they have been widely studied at hadron 
colliders. Anomalous couplings have been studied in top pair production izHU!, top pair plus jets [15], direct photon 
production m and single top production m- Also, spin correlations in top pair production have been used in the 
study of CMDM[T5|. Constraints to CMDM and CEDM from Higgs boson production at the LCH have also been 
reported CD- 

Anomalous moments have been calculated in different new physics scenarios as in the case of MSSM |20| , THDM [21] , 
Little Higgs model [22] and Unparticles [23]. CMDM and CEDM are dehned through the effective Lagrangian 

C = (Ak + iqsAd) u{t), (1) 

where Ak and Ad represent the CMDM and CEDM respectively, is the gluon field strength, and T“ are the 
QCD fundamental generators of SU(3)c. 
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In a recent study by the CMS collaboration spin correlation in ti cross section is used to obtain a new bound to 
CMDM; the bound found is Re(A«;) = 0.037 ± 0.041, at 95%CL, or equivalently —0.045 < Re(AK;) < 0.119 |24) . 
The given bound is obtained in the dilepton channel in pp collisions at y/s = 7 TeV with an integrated luminosity of 
5 fb“^ and is compared with the SM theoretical prediction, including a new physics contribution. Bounds to the top 
quark CMDM and CEDM are found from combined results from Tevatron and LHC (Atlas) using high values of the 
rritt cross section, the bounds reported are |Ak| < 0.05 and Ad < 0.16 at 95%C.L. [B]. It is estimated that in order 
to find a 5 (t upper bound of Ad < 0.05 it is required an integrated luminosity of 10 fb“^ at \/s = 13 TeV at LHC [5]. 

From CLEO data in the radiative decay 6 —>■ sq a strong constraint of —0.03 < |Ak| < 0.01 is obtained [^. The 
SM prediction to the CMDM is Ak ~ 5.6 x 10“^ [2S]. The top quark CEDM and CMDM induce new contributions to 
the lightest quarks through the renormalisation group equations (RGE), therefore, the neutron dipole moment gives 
an indirect bound to CEDM |Ad| < 1.9 x 10“^ [ 6 ]. 

In the present work we want to study the THDM Type-HI contribution, without imposing a Z 2 discrete symmetry, 
in the quark top CMDM and CEDM. This kind of model explicitly violates the CP symmetry in the scalar potential, 
which generates the top quark CMDM at one loop level. Explicit CP violation generates mixing among the neutral 
CP-odd and CP-even scalar fields. This mixing is very suppressed by the recent data obtained for in LHC, for 
the Higgs physics [5B]. Using the allowed region ai — a 2 from the neutral scalar sector we find the allowed values for 
the top quark CDMD and CEDM. 


II. TWO HIGGS DOUBLET MODEL WITH CP VIOLATION 


The simplest extension of the SM, with one extra scalar doublet is called the two Higgs doublet model (THDM); the 
model contains two doublet fields, <i)i and $ 2 ; this kind of model has the advantage of being capable to describe the 
phenomenon of CP violation [28]. When a discrete symmetry is imposed there are two kind of models, in the so called 
Type-I one doublet gives mass to all quarks and in the Type-H model one doublet gives mass to the up quarks while 
the other one gives mass to the down quarks. In a theory without the restriction of a discrete symmetry, also called 
THDM Type HI, the two doublets simultaneously give mass to the up and down quarks, and the mass matrix depends 
on the Yukawa couplings which cannot be simultaneously diagonalized, allowing the presence of flavour changing at 
tree level |^ . 

If we consider a general THDM the scalar potential can be written as 


V 


$1 - $2 - + h.c] 

^Ai -b ^A 2 ($^$ 2 )% A 3 (4>^4>i) ($^$ 2 ) + A 4 (4>^4’2) 

2^5 (‘^’^^ 2 ) -b Ag (4)^$i) (<I>^$ 2 ) -b Ay ($^<I> 2 ) (4)]^$2) -b h.c. , 


( 2 ) 


where pi, /i|, Ai, A 2 , A 3 and A 4 are real parameters and the parameters •^ 5 j ^61 and Ay can have complex values 
allowing the explicit CP violation in the potential. The neutral components in the fields are dehned as -^{va+rja+iXa): 
where a = 1,2. The vacuum expectation values (VEV) can be taken real because complex phases can be reabsorbed 
in the complex parameters in the scalar potential. The VEV take the values 


(‘fi) 


1 



(3) 


($2) 


1 

71 



(4) 


Due to the explicit CP symmetry breaking, there will be mixing among the CP-odd and CP-even scalar sectors. 
Defining tan/3 = we take the scalar field ( 7/3 = —XiS/s -b orthogonal to the Would-be Goldstone component 
corresponding to tfie Z gauge boson. After symmetry breaking, the mass eigenstates of the neutral Higgs bosons are 
related to the rjj states as 


3 

h j — ^ ^ Vj : 


(5) 
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where i = 1, 2, 3 and the R matrix is given by |29j : 

( ClC 2 SiC 2 32 \ 

- (C1S2S3 + S1C3) C1C3 - S1S2S3 C2S3 , (6) 

-C 1 S 2 C 3 + S 1 S 3 - (C 1 S 3 + S 1 S 2 C 3 ) C 2 C 3 J 

where the abbreviations Ci = cosai and Si = sina^, with i = 1,2,3 are used, hi eigenstates do not have a well defined 
CP state. For convenience, we choose the neutral Higgs bosons hi to satisfy the mass hierarchy irihi ^ ^ fnhi- 

In the limit case when S 2 = S 3 = 0 we recover the THDM without CP violation. 

The Yukawa Lagrangian for the quark sector has the general form 

3 2 

- Cvuka^a = E E + h.C.) . ( 7 ) 

iij=l a—1 

In the above equation, are the 3x3 Yukawa matrices, ql denotes the left handed quark doublets and ur: 

da, represent the right handed quark singlets under SU{2)l. After spontaneous symmetry breaking (SSB), the mass 
matrix can be written as 


2 

^ 72 “ 

n=\ ’ 


( 8 ) 


where YJ = V[Y^h rotation matrices that diagonalize the mass matrix. The 

Yukawa matrix Y 2 as a function of M“ and Y“ gives THDM-II Lagrangian, with tree level flavor changing. For the 
up sector the Yukawa Lagrangian can be written as 

ijk 

p .., 

Ilk 

(9) 

where 

= Rk 2 -iRks cos 13, 

Bl = i?fcisin/3-i?fc2COS,0 + ii?fc3. (10) 


The Yij also gives a contribution to the anomalous couplings CEDM and CMDM wich is of the same order of the 
one in the THDM-II coming from {A^Pr -|- A^Pr) Ujhk). If we asume the Cheng-Sher parametrization |30j . 

where (Yu ~ rrit/v), both contributions must be taken into account in order to compute Ak y Ad. 


III. CMDM AND CEDM IN THE GENERAL THDM 

The anomalous couplings contributions for the CMDM Ak and for the CEDM delta Ad, arising from the diagram 
in Fig. 0 are given by the following expressions 


Ak = 




3 p\ 

E / 


dy-. 


2727r2 sin^ Jo {x + y)'^ - [x + y - l)™^ 

[(x -h y){x + y- l)(i?^2 - cos^ ^Rf^) - {x + y){Ra + cos^ JRa)] > 


( 11 ) 


Ad = - 


Gpmf 


sin^ /3 
cos jdRaRa], 


3 

E 

2=1 ’ 


/•i 

dx dy 

0 Jo 


{x + y){x + y-l) 

{x + y)"^ - {x + y- 1)to? 


( 12 ) 
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FIG. 1: Feynman Diagram for the anomalous quark-gluon couplings in the general THDM. 


where mi, m 2 , and m 3 are the masses of the hi, h 2 , and h^, respectively. In this calculations we have used the method 
presented in [illlllllS]. 


The contribution from the Yu is 


Akm = 


Gpru^ 


3 pl — x 


H / dx 


dy-, 


1 


2^2712 sin^ Pj^iJo Jo {x + y)^ -{x + y- l)ml^ 


(13) 


{x + y){x + y - 1) {{Rii sin/3 - cos/3)^ - Rj^) - {x + y) {Rn sin,5 - Ra cospf + R% 

and for the CEDM we have 

Gjrm? ^ 


Adtt — — 


^ / (ii 


<is—X 


y/2TT^ sin^ P Jo {x + yy -{x + y- l)m; 


(14) 


( Rii sin f3 - Ri 2 cos (3)Ri 3 . 


The contributions for CMDM and CEDM from the coupling proportional to M“ in one vertex and F“ in the other 
vertex are given by 


A^^ni — 


Gprul 


3 „1 




f'l — X 


dy-. 


1 


2v^7r2 sin^ /3 ^ ^0 Jo (x + y)'^ - (x + y - l)ml^ 
{x + y){x + y - 1) ((i?ii sin/3 - Ri 2 cos/ 3 )i?i 2 + Riscos/S) 
-{x + y) {{R,i sin /3 - R^ cos /3)Ri2 - R% cos ji) 


(15) 


and 


^dir}.t — 




H / dx 


V2Tr^ sin^ /3 ^ 7o Jo {x + y)^ - {x + y - l)m; 


(16) 


{x + y){x + y - 1) {RaRa - {Rn sin/3 - Ra cos j3)Ra cos/3) 
-{x + y) {Ri 2 Rto + {Rii sin/3 - Ri 2 cos,5)13*3 cos(3) 


We are using the Cheng-Sher parametrization of Yu = mtiv in Eqs. (11) - (16). We denote three contributions as 
Akj = Ak -|- Aiftt -I- Akint and Adt = Ad -I- Adu + Admt- The charged Higgs contribution to Ak and Ad can be 
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neglected because in the loop circulates a bottom quark and it is suppressed compared to the loop contribution in the 
neutral scalar sector where a top quark is the circulating particle. For the second and third terms in the sum of Eqs. 
@ @ we considered mh^ = = ^h+ ■ In this partially degenerate case there is CP violation for a 2 ^ 0 m- 

We will study nine regions of interest in the ai-a 2 parameter space, this approximate regions are described in Table 
|l] and are those already under consideration by the authors in previous work |32j . The allowed regions Ri..g in the 
ai — 02 plane, are defined from experimental bounds in |27j . wherei?-^.^ is given by 


o'(gg ^ hi)Br(hi -» 77) 

a(gg hsM)Br{hsM 77)' 


(17) 


Because the charged Higgs contributes to the loop in hi 77 , the chosen values of Mh± and tan/3 affect the allowed 
region in ai — 02 - The process {B —>■ Xs"/) contains an important contribution from the charged Higgs, this process 
strongly restricts Mu± vs. tan/3 [55]. For small values of tan/3 the bound to the charged Higgs mass is of around 
300 GeV [321 ■ A global analysis of B decays restricts Mjj± < 400 GeV and tan/3 < 10 |34H55] . 

In Table |l| are shown the Ri regions for the given values of M^+ and tan /3. In each region we set the masses of the 
neutral Higgses rrih^ and rrih^ equal to the mass of the charged Higgs = nT-ha = . 

Using (0.5 < R^y < 2.0), Mjj± = 300 GeV and tan/3 = 1.0, the allowed regions in the ai — 02 plane are i?i and 
i? 2 . For the same values in the other parameters and setting the charged Higgs mass to = 500 GeV region 

i ?3 is obtained. Combining (1 < Rry^ < 2.0) with = 350 GeV and tan/3 = 1.5 the allowed region for ai — 02 

give i ?4 and R^. When tan/3 = 2 and for the same values to the other parameters as in i ?4 5 we get i? 6 , 7 . Finally if 
tan/3 = 2.5 the allowed regions are i? 8 , 9 . 

In Table we present the limits obtained for Xkt and Ad* (or Ak and Ad) in a general THDM. In order to 
illustrate the limits reported in Table |h] we show in Figs. (§, (§ and Q the limits obtained for the anomalous 
moments, with values of the alpha parameters allowed for each region. In order to estimate the contribution of each 
term we separately analyze different cases. In the Fig. ( |^ o nly the contribution from Eqs. (11) and is considered, 
meanwhile in Fig. ([^ the contribution from V“, Eqs. (|13| ) and (14), to previous values shown in Fig. ([^ is added. 
All the contributions, including the interference terms (15) and (16), are considered in Fig. 
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FIG. 2: Scatter plot in the and Ad plane with random values of the angles ai, 0:2 in the range allowed for each region and 
0(3 = 0., ta.nP and are as defined in Table Q in every region. 
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FIG. 3: Scatter plot in the Akt and Adt plane with random values of the angles Qi,Q 2 in the range allowed for each region 
and as = 0., tan/3 and Mu+ are as defined in Table (j^ in every region. In this case we plot Adt and ARt including the Yu 
contribution. 
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FIG. 4: Scatter plot in the ARt and Adt plane with random values of the angles Qi,Q 2 in the range allowed for each region 
and as = 0., tan/3 and Mu+ are as defined in Table (j^ in every region. In this case we plot Adt and ARt including the Yu 
contribution and also the interference contribution. 
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TABLE I: M^f+and tan/3 in each region. 



«1 

02 

M^+(GeV) 

tan j3 

Ri 

0.67 < oi < 0.8 

0 < 02 < 0.23 

300 

1 

R2 

0.8 < ai < 1.14 

-0.25 < 02 < 0. 

300 

1 

Rs 

1.18 < ai < 1.55 

-0.51 < 02 < 0 

500 

1 

R 4 

-1.57 < ai < -1.3 

-0.46 < 02 < 0. 

350 

1.5 

R 5 

0.93 < ai < 1.57 

-0.61 < 02 < 0. 

350 

1.5 

Re 

-1.57 < ai < -1.28 

-0.38 < 02 < 0. 

350 

2 

R 7 

1.08 < ai < 1.57 

-0.46 < 02 < 0. 

350 

2 

Rs 

-1.39 < ai < -1.3 

-0.13 < 02 < 0. 

350 

2.5 

Rg 

1.16 < ai < 1.5 

-0.43 < 02 < -0.1 

350 

2.5 


TABLE II: Range of values taken by the anomalous CMDM and CEDM. The first row shows only the contributions of Ait 
and Ad according to Fig. The regions in Fig. |^, which correspond to the addition of the contribution are shown in 
second row. The last row shows the regions for the three contributions based in Fig. 


R 

CMDM 

CEDM 

Ri 

-1.42 X 10“^ <Ak< -1.33 X 10“^ 

0.< Ad< 2.41 X 10""^ 


-2.73 X 10~^ < Akt < -2.63 x 10"^ 

0. < Adt < 4.56 X lO"'^ 


-2.76 X 10~^ < Akt < -2.50 x 10"^ 

-1.26 X 10“® < Adt < 0. 

R 2 

-1.67 X 10“^ <Ak< -1.42 x 10“^ 

-5.30 X 10“'‘ < Ad < 0. 


-3.06 X 10~^ < Akt < -2.37 x 10"^ 

-5.15 X 10“'‘ < Adt < 0. 


-2.50 X 10“^ < Akt < -1.67 x 10"^ 

0. < Adt < 2.65 X 10"^ 

Re 

-1.72 X 10"^ <Ak< -1.40 X 10“^ 

-1.61 X 10“® < Ad < 0. 


-3.03 X 10~^ < Akt < -2.43 x lO"^ 

-6.63 X lO-'^ < Adt < 0. 


-1.35 X 10“^ < Akt < -6.25 x 10"® 

0. < Adt < 5.13 X 10"^ 

R 4 

-1.18 X 10~^ <Ak< -1.11 X 10“^ 

0. < Ad < 7.48 X lO""^ 


-2.35 X 10“^ < Akt < -2.14 x 10"^ 

-3.12 X 10“'‘ < Adt < 0. 


-9.84 X 10“® < Akt < -5.91 x 10"® 

-4.36 X 10“® < Adt < 0. 

Re 

-8.89 X 10“® <Ak< -1.22 x 10“^ 

-7.36 X 10“'‘ < Ad < 0. 


-1.79 X 10“^ < Akt < -2.23 x 10"^ 

-8.23 X 10“'‘ < Adt < 0. 


-1.80 X 10~^ < Akt < -6.33 x 10“® 

0. < Adt < 5.15 X 10"^ 

Re 

-9.83 X 10“® <Ak< -9.59 x 10“® 

0 < Ad < 4.49 X lO”'^ 


-1.94 X 10'^ < Akt < -1.79 x 10"^ 

-2.69 X 10“'‘ < Adt < 0. 


-9.71 X 10“® < Akt < -5.62 x 10“® 

-3.45 X 10“^ < Adt < 0. 

Rr 

-1.03 X 10'^ <Ak< -8.30 x 10“^ 

0 < Ad < -4.78 X lO"'^ 


-1.84 X 10~^ < Akt < -1.58 x 10"^ 

-5.11 X lO-'^ < Adt < 0. 


-1.59 X 10“^ < Akt < -6.99 x 10"® 

0. < Adt < 3.98 X 10"^ 

Rs 

-9.13 X 10“® <Ak< -9.01 X 10“® 

0. < Ad < 1.30 X lO""^ 


-1.74 X 10“^ < Akt < -1.71 X 10"^ 

-9.02 X 10“® < Adt < 0. 


-7.93 X 10“® < Akt < -6.55 x 10"® 

-1.23 X 10“® < Adt < 0. 

Rg 

-9.30 X 10~® <Ak< -7.87 x 10“® 

-3.65 X lO-'^ < Ad < -1.04 x 10“'‘ 


-1.63 X 10“^ < Akt < -1.48 x 10"^ 

-3.97 X 10“'‘ < Adt < -1.84 x 10“® 


-1.44 X 10'^ < Akt < -8.0 X 10“® 

8.31 X lO"'^ < Adt < 3.59 x 10"® 


IV. CONCLUSIONS 

In this work we have studied regions of interest in the ai —a^ parameter space, in order to calculate the contribution 
to the top anomalous couplings CMDM and CEDM in the context of a general THDM with CP violation. In 
our analysis has been considered the contribution of the Yukawa coupling Yu using M//+ = 300,350 GeV and 
tan/3 = 1,1.5,2,2.5. We find for the nine regions of interest that the value for Ak can be at most of order 10“^ and 





























Ad of order 10“^ . The contributions arising from the interference of M“ and Yu have been considered in the results. 
The contributions of Akj are added coherently and for the three contributions the variations are not appreciable for 
the different regions Ri, i = 1, However the contribution of Adt coming from the interference increases almost 
an order of magnitude for the regions Ri and in some cases the sign is changed. 

The recent NLO calculation for top-quark production including anomalous top-quark CMDM reports the bound 
—0.0096 < Ak < 0.0090 [13]; our result for region Rr is in agreement with this stringent constraint. We can 
also compare our theoretical bounds with those obtained from Higgs boson production at the LHC where more 
conservative model independent bounds are obtained m. our results in all nine regions are in agreement even with 
the most restrictive bounds projected at 14 TeV, —0.016 < Ak < 0.008 and |Ad| < 0.007, as reported in Table [h| of 
reference m- 

With future LHC measurements at higher energy there will be an excellent chance to probe new physics properties 
of the top quark; anomalous dipole moments are a measure of this new physics properties that can also give some 
insight in the top quark structure. A precise measurement of the top quark CMDM and CEDM, expected soon after 
future LHC results, will be a useful source of information in order to discriminate among different SM extensions. 
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